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H ' Abstract 

o: 

, We study a general class of nonlinear second-order variational inequalities with interconnected bilateral 

obstacles, related to a multiple modes switching game. Under rather weak assumptions, using systems of 
penalized unilateral backward SDEs, we construct a continuous viscosity solution of polynomial growth. 
I^H ■ Moreover, we establish a comparison result which in turn yields uniqueness of the solution. 
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1 Introduction 

In this paper we study systems of variational inequalities with interconnected lower and upper obstacles. This 
type of inequalities arises as the Bellman- Isaacs equation in a multiple modes switching game between two players. 
Besides their classical fields of applications, multiple modes switching games are attracting a lot of interest in the 
0^ ■ management of power plants (see Bernhart (2011) ([3) and Perninge (2011) ([27]), where they are successfully 
used to design optimal stopping and starting strategies for power flow control through activation of regulating 



bids on a regulated power market. 



£SJ \ The objective of this work is to establish existence and uniqueness of a continuous viscosity solution of the 
following system of variational inequalities with oblique reflection: 



mm{(vV -L ij [v\)(t,x), max{(^' -U ij [v\)(t,x), 

-d t v^(t, x) - (i, x) - /«(t, x, (v kl (t, a:)) (fc , I)e ri x r 2 , <r T (i, x)D x v^(t, x))} } = 0, (1.1) 
v*i(T,x) = h ij (x), 

for every pair in the finite set of modes T 1 x V 2 , where, for any (t,x) £ [0, T] x M fe , 

&p(t,x) := b(t,x)D x cp{t,x) + -Tr[a<r T (t,x)Dl x ip(t ) x)], 
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and to any solution v — (u t5 )(t,j)er 1 xr a we associate the lower obstacle operator 

L ij [v\(t,x) := max (v kj - g..)(t,x), 

and the upper obstacle operator 

U^M(t,x) := min 

where, g (resp. ^7) stands for the switching cost incurred when the first (resp. second) player decides to switch 
from mode i to mode k (resp. from mode j to mode I). Finally, the function f lJ stands for the instantaneous 
payoff when the first player is in mode i and the second one in mode j. 

The system (jl.ip and related switching games have been studied by several authors. The most recent work 
discussing this topic includes the papers by Hu and Tang (2008) ([H]) and Tang and Hou (2007) ([IE]) ( see a l so 
the references therein) which deal with switching games related to p.l[) . when the switching costs do not depend 
on the state variable. To the best of our knowledge, Ishii and Koike (1991) ([H]) are the latest most general 
existence and uniqueness results for the system (|1.1[) , for state-dependent switching costs. They derive existence 
and uniqueness of viscosity solutions of the elliptic version of (|1.1[) in a bounded domain of M fe whose boundary 
is of class C 2 , when the so-called Fichera functions are strictly negative (see [18 , Proposition 4.3). 

The main result of the present paper, which is given in Theorem 15.41 establishes existence and uniqueness of a 
continuous viscosity solution of the system (|1.1[) . when the state space is the whole M. k and under rather weak 
assumptions on the involved coefficients. Our approach is probabilistic and makes use of penalization schemes 
that allow us to connect the related penalized PDEs with systems of reflected backward SDEs with unilateral 
interconnected obstacles which, for instance, have been studied in [7], [11], [17] or [20]. With the help of these 
sequences of solutions of reflected BSDEs and their connection with PDEs, via Feynman-Kac's formula, we are 
able to construct in Propositions 15.11 and 15.21 both a viscosity subsolution and a supersolution for the system 
(jl.ip . Relying next both on the comparison result established in Theorem 14.21 and adapting the Perron's method 
we construct a solution for (jl.ip which is therefore unique. Finally, using again the uniqueness result, we identify 
the limit of the penalized decreasing scheme as the solution of the system (jl.ip . 

We made this detour instead of trying to solve a related system of reflected BSDEs with interconnected bilateral 
obstacles, as one would expect, simply because they satisfy neither the so-called Mokobodski condition nor 
the condition of strict separation of the two obstacles, which would guarantee existence and uniqueness of the 
solutions of the system of BSDEs, since these obstacles depend on the solution. The structure of these bilateral 
obstacles suggests a rather new type of conditions to guarantee existence and uniqueness result for the related 
system of reflected BSDEs. This problem is beyond the scope of the present paper and therefore left for future 
research. 

Our plan for this paper is as follows. In Section 2, we provide all the notations used in the paper, state the whole 
list of required assumptions and define viscosity sub- and supersolutions along with equivalent characterizations. 
In Section 3 we construct two approximation schemes (an increasing and a decreasing on), consisting of sequences 
of penalized reflected BSDEs associated with standard switching problems. The counterpart of the decreasing 
scheme (resp. the increasing one) in terms of PDEs stands for the penalized scheme of system (jl.ip (resp. system 
(jl.2p given below). Section 4 is devoted to the proof of a comparison result related to the sub- and supersolutions 
of (jl.ip . In Section 5, the decreasing limit is identified as a viscosity subsolution of (jl.ip while a super-solution is 
exhibited. Finally, we use Perron's method to construct a viscosity solution of (jl.ip and, thanks to the uniqueness 
result, its connection with the limit of the decreasing scheme is obtained. As a by product, we also obtain a 
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similar characterization of the limit of the increasing scheme as the unique solution in viscosity sense of the 
following system of variational inequalities: For every £ T 1 x T 2 , 

{max { (v 1 ? - U ij [v\) (t, x) , min { (v ij (t, x) ~ \v\) (t, x) , 
-dtvv (t, x) - Cv 13 (t,x)- P (t, x, (v kl (t, x)) ( fc,i )e ri xP,^ (t, x)D x v» {t,x))}}=0, (1.2) 
v ij (T,x) = h ij (x). 

We do not know whether the solutions of (jl.ip and (ll.2j) coincide or not. We note that this is a very important 
issue since this will enable to characterize this common solution as the value of the zero-sum switching game. 
Since it is beyond the scope of the paper, this question is left for further research. 

2 Preliminaries and notation 

Let T (resp. k, d) be a fixed positive constant (resp. two integers) and T 1 (resp. T 2 ) denote the set of switching 
modes for player 1 (resp. 2). For later use, we shall denote by A the cardinal of the product set T 1 x T 2 and for 
G T 1 x T 2 , (r 1 )" 1 := r 1 - {z} and (r 2 p := r 2 - {j}. Next, for y = (y w ) (fc)I)er i xra G M A , G r 1 x r 2 , 
and y G R, we denote by \y~^\ y] the matrix which is obtained from y by replacing the element y % i with y. 

Next, let us introduce the following functions. For any G T 1 x T 2 , 

b : (t, x) G [0, T] x M l 4 b{t, x) G K fe ; 
a : (t, x) G [0, T] x R fc h-> a(t, x) G R kxd ; 
P : (t, x, y, z) G [0, T] x R fc+A+d ^ (i, x, y, z) G K ; 
(t,x)e[0,T}xR k ^gJt,x)eR (k G (r 1 ) - *); 
: (i,s) G [0,T] xR^ g 3 i{t,x) el (Z G (r 2 )-J); 
h ij : x G R fc i-> /^(x) G K. 

A function $ : (i, a;) G [0, T] x M 70 H- $(i, 2) € R is called of polynomial growth if there exist two non-negative 
real constants C and 7 such that 

|$(i,x)| < C*(l + |x| 7 ), (t,x) G [0,T] x M fe . 

Hereafter, this class of functions is denoted by Let C 1,2 ([0,T] x R fc ) (or simply C 1 ' 2 ) denote the set of real- 
valued functions defined on [0,T] x M fc , which are once (resp. twice) differentiable w.r.t. t (resp. x) and with 
continuous derivatives. 

In this paper, we investigate existence and uniqueness of viscosity solutions v(t, x) := (v kl (t, x))( k ^ er i xF 2 of the 
following system of variational inequalities with upper and lower interconnected obstacles: For any G T 1 x T 2 , 

min { (v ij - L ij [it]) (t, x), max { (v^ ~ U lj [v\) (t, x), 

-dtv* (t, x) - Cv* (t, x) - P (t, x, (v kl (t, x)) (M)£rl xF 2,a T (t, x)D x v^ (t,x))}}=0 (2.1) 
v ij (T,x) = h ij {x) 

where, for any (t,x) G [0,T] x M fe , 

£<p(t,x) := b(t,x)D x <p(t,x) + ^Tr[aa T (t, x)D 2 xx ip{t, x)), 

and G T 1 x T 2 , 

L ij [v\(t,x) := ^ max_ {v k3 {t,x) -g ik (t,x)) and U %3 [v\(t,x) = f min_ (v il (t, x) + g jt {t, x)). 
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The functions f lJ stand for the instantaneous payoff when the first player is in mode i and the second one in 
mode j, and g (resp. g~ji) stands for the switching cost incurred when the first (resp. second) player decides to 
switch from mode i to mode k (resp. from mode j to mode I). 

The lower obstacle L %] [v\ and an upper obstacle U 1 ^ [v] are called interconnected because each of them depends 
on the underlying solution v :— (t |fe ')(fc,i)sr 1 xr 2 - 

In a way, the system (I2.1[) is the Bellman-Isaacs system of equations associated with the zero-sum switching game 
with utility functions (/ y )(i.j)er 1 xr 2 i terminal payoffs (/& ,J )(t j^er 1 xr 2 and switching costs for the maximizer 
(resp. minimizer) given by (g..) {iJ)er i xr 2) (resp. (<7iy)(j ,j) e r> xr 2 ))- 

The following assumptions are in force throughout the rest of the paper. 

(HO) The functions b and a associated with the second order operator C are jointly continuous in (t, x), of linear 
growth in (t, x) and Lipschitz continuous w.r.t. x, meaning that there exists a non-negative constant C 
such that for any (t, x, x') G [0, T] x R k+k we have: 

|&(t,x)| + \<r{t,x)\ < C(l + \x\) and \a(t,x) - a(t,x')\ + \b(t,x) - b(t,x')\ < C\x - x'\. 
(HI) Each function P 

(i) is continuous in (t, x) uniformly w.r.t. the other variables (y, z) and for any (t, x) and the mapping 
(t, x) — > f l i (t, x, 0, 0) is of polynomial growth. 

(ii) satisfies the standard hypothesis of Lipschitz continuity with respect to the variables (y := (y u )(ij)eri xr 2 , z ), 
i.e. V (t,x) £ [0,T] x R fc , V (yi,y 2 ) £ M A x R A ,(z\z 2 ) £ R d x R d , 

- f j (t,x,y 2 ,z 2 )\ <C(\y x -yi\ + \z x -z 2 \), 

where, \y\ stands for the standard Euclidean norm of y in R A . 

(H2) Monotonicity: Let y = (y kl )(k,i)er^xr 2 , then for any e L 1 x L 2 and any (k,l) ^ (i,j) the mapping 

yk,l _^ f l i(s,y,z) is non-decreasing. 

(H3) The functions h %0 , which are the terminal conditions in the system (|2.ip . are continuous with respect to x, 
belong to class II g and satisfy 

V (i,j) £ T 1 x L 2 and z e R fc , ^ max . (/i fej (a;) - g. k (T,x)) < h l3 (x) < ie ^_ 3 (h u (x) + gji(T,x)). 

(H4) The no free loop property: The switching costs g ik and g~ji are non- negative, jointly continuous in (t, x), 
belong to H g and satisfy the following condition: 

For any loop in T 1 xF 2 , i.e., any sequence of pairs (i\, j\ ),..., (in, Jn) of T 1 xT 2 such that (In, 3n) — (h,ji), 
card{(ii, ji), . . . , (ijv,iiv)} = N — 1 and \/ q = 1,...,N — 1, either = i q or = j g , we have 

V(t,x) G [0,T] x R fc , 

q=l,N-l 

where, V g = 1, . . . , N - 1, ^ i (t, x) = -g (t, x)l iq7 t iq+1 + g jqig+1 (t, x)t jq ^ q+1 

~"q fc q + 1 

(resp. ip lqlq+1 (t,x) =g 4 ^(t, x)t lq ^ lq+1 - g jq , lq+1 (t, x)\ Jq -L jq+1 ). 
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This assumption implies in particular that 

N-l 

V (ii, In) G (r 1 )^ such that ijv = ii and card{ii, . . . , In} = N — 1, 75. . >0 (2.3) 

p=l 

and 

N-l 

V (ji, • • ■ Jn) 6 (r 2 )^ such that jjv = ii and card{ji, . . . , j N } = N - 1, ^ 9j h ,j k+1 > 0. (2.4) 

P =i 

By convention we set g~jj = g . . = 0. 

Conditions (|2.3I) and (|2.4p are classical in the literature of switching problems and usually referred to as the no 
free loop property. 

Finally, let us mention that if we set 

9 u (t,x) = \i-j\g(t,x) and g~ij(t,x) = \i - j\g(t,x), (i, j) G T 1 x T 2 , 

where both g and g are functions such that, for any (t,x) G [0, T] x M fc , is not a rational number, then 

assumption f|2.2[) holds. 

We now define the notions of viscosity super (or sub)-solution of the system (|2.1j) . This is done in terms of the 
notions of subjet and superjet which we recall here. 

Definition 1. (Subjet and superjet) 

(i) For a lower semicontinuous (Isc) (resp. upper semicontinuous (use)) function u : [0,T ] x R k ->■ E, we denote 

by J~u(t,x) (resp. J + u(t,x)) the parabolic subjet (resp. superjet) of u at (t,x) G [0, T] x E fe , as the set of triples 
(p, q, M ) G E x M fe x § fc satisfying 

u{t',x') > {resp. <) u(t,x) + pit! -t) + q T (x' - x) + 

\(x' - xYM{x' -x) + o(\t' -t\ + \x' - x\ 2 ) 

where S fe is the set of symmetric real matrices of dimension k. 

(ii) For a Isc (rep. use) function u : [0,T] x R fc — > E, we denote by J~u(t,x) (resp. J + u(t,x)) the parabolic 
limiting subjet (resp. superjet) of u at (t,x) G [0, T] x WL k , as the set of triples (p,q,M) G E x R fe x § fc such that: 

(p,q,M) = lim n (p„,q„,M„), (t,x) = lim n (i„, x n ) with [p n ,q n ,M n ) G J~u(t n ,x n ) 
(resp. J + u(t n ,x n )) and u(t, x) = lim n u(t n , x n ). 



Finally, given a locally bounded K-valued deterministic function u, we denote by it* (resp. u*) its lower (resp. 
upper) semicontinuous envelope defined by: 

V (i, x) G [0, T] x K fc , u*(t,a;)= lim u(i',x') and ii*((,i)= lim u{t',x'). (2.5) 

(t' ,x' )^(t,a:); t'<T (t',x')->-(t,3;); t'<T 

We now give the definition of a viscosity solution for the system (|2.ip . 



Definition 2. Viscosity solution to H2.1\) 

(i) A function v — (v kl (t, x))^ k ^ er i xT 2 ■ [0,T] x R k n> M A such that for any (i,j) G T 1 x T 2 , v lj is Isc (resp. 
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use), is called a viscosity super solution (resp. a viscosity subsolution) to \2.1\) if for any (i,j) G L 1 x V 2 , for any 
(t,x) G [0,T) x R fe and any (p,q,M) G J~v tJ (t,x) (resp. J + v l3 (t,x)) we have: 

min x) — L tJ '[v\(t, x), 

max{-p-b(t,x).q - \Tr[{aa T )(t,x)M] - fi (t, x, v(t, x), cr T (t, x)q); v ij (t, x) - U^[v\(t,x)}} > 0, 
Vi(T,x) > h^(x), 

(2.6) 

(resp. 

min {v lJ (t, x) — [v\(t, x), 

max.{-p-b(t,x).q - ±Tr[(acr T )(t,x)M} - f^(t,x,v(t,x),a T (t,x)q);v ij (t,x) - U lJ [v\(t,x)}} < 0, 
v ij (T,x) < h ij (x)y 

(2.7) 

(ii) A locally bounded function v = (t |fe ')(fc,/)er 1 xr 2 : [0,T] xR 1 -) M A is called a viscosity solution of \2. 1}) if the 
associated lower (resp. upper) semicontinuous envelope (v* 1 )^ i)er 1 xr 2 ( res P- i vkl )(k,l)ev 1 xr 2 ) defined in 12.5]) 
is a viscosity supersolution (resp. subsolution) of 12. 1\) . 

If, in addition, for any (k,l) G L 1 x T 2 , v kl — v M * , then (v kl )^.i)er 1 xr 2 * s a continuous viscosity solution of 

Remark 1. Under Assumptions (H0)-(H4), the above definition of a viscosity solution for \2.1)) can be relaxed 
replacing the limiting subjet J~(vl J )(t,x) of the supersolution v* (resp. the limiting superjet J + (v lJ '*)(t,x) of v* ) 
by the subjet J~(vl)(t, x) (resp. by the superjet J + (v l ' h *)(t, x)). This results from the continuity of the functions 
b, a, (f 3 ,h' i: ',g. j ,gij)^ j ) er i x r^ and the monotonicity property (H2) o/ (/^) (iJ)er i )<r 2. 

3 Systems of reflected BSDEs and approximation schemes of the 
solutions 

Let (f2, T, P) be a fixed probability space on which is defined a standard d-dimcnsional Brownian motion B = 
(B t )o<t<T whose natural filtration is (.F t := cr{B s , s < i})o<t<T- Let F = (J r t)o<t<T be the completed filtration 
of (J-f )o<t<T with the P-null sets of J 7 , hence (J-"t)o<t<T satisfies the usual conditions, i.e., it is right continuous 
and complete. Furthermore, let 

• V be the cr-algebra on [0, T] x of F-progressively measurable sets; 

• H 2 ' e (£ > 1) be the set of "P-measurable and Revalued processes w = (wt)t<T such that ^[Jq \w s \ 2 ds] < oo; 

• S 2,i (£ > 1) be the subset of H 2 e of continuous processes such that E[sup t<T \w t \ 2 ] < oo. Finally let A + ' 2 
be the subset of S 2 ' 1 of non-decreasing processes K = (K t )t<T such that K = and E[if|,] < oo. 

Next, for n, m > 0, let (Y ij ' n ' m , Z^"< m ) (lj)er i x r 2 be the solution of the following system of BSDEs. 

dY i3,n..m = _/^«^( S; X*' K ,(y s fe/ '"' m ) (fc , 0er i xr 2,Z^ n ' m )ds + Z^"' m dS s , S<T, (3.1) 

Y*> nm = h i i{x% x ), 

where, 

fij,n,m^ x t,^ {y i j)( . j)erixr2;Z) ■.= fi( s ,X t s '-,(y k % il)erlxr2 ,z)+ 

<V lJ - max fe(rl) -,{^ - g ik (s, Xp*)})- - m(y ij - mm le(r2) -i {y il + g jt (s, X s *- X )})+. 
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Note that by Assumption (HI) and the standard result on multi-dimensional BSDEs by Pardoux and Peng (1990) 
(ED), the solution (y y, ™' m , Z lJ ' n ' m )^^ er i xr 2 exists and is unique. Since the generators and the terminal values 
depend on (t,x), the processes Y lJ ' n ' m and Z l3 > n ' m also depend on (t,x) but, to avoid overload notation, we do 
not mention this dependence in the sequel. Furthermore, the following monotonicity properties holds for the 
double sequence (Y lJ > n > m ) nym , 

Proposition 3.1. For any G T 1 x T 2 and n,m > we have 

P-a.s., Y ij < n ' m < Y ij > n+1 ' m and Y ij ' n ' m > Y ij > n > m+1 , (i, j) G T 1 x T 2 . (3.2) 

Moreover, for any (i,j) G T 1 X T 2 and n, m > 0, there exists a deterministic continuous function v l °' n,m in H g 
such that, for any t < T, 

Y»' n ' m =v ii ' n ' m (s,X t s ' x ), se[t,T]. (3.3) 

Finally, for any (i,j) £ T 1 X T 2 and n,m> 0, 

v ij ^ m (t,x) < v^ n+1 - m (i,x) and v ij ' n ' m (t,x) > v l3 ^ m+1 (t, x), (t,x) G [0,T] x M fe . (3.4) 

Proof. The second claim is just the representation of solutions of standard BSDEs by deterministic functions 
in the Markovian framework (see e.g. El Karoui et al. (1997) ([I3j) for more details). As for the first one, it 
is based on the result by Hu and Peng (2006) ([22]) related to the comparison of solutions of multi-dimensional 
BSDEs (we recall in Appendix (Al)). To prove this, it is enough to show that for any t, (yij)(%,j)£r 1 xr 2 -i 
(yij)(<j)ePxP S K A and (zij) (i j )er i xF 2, (%) (ij)er i xr2 6 (K d ) A we have: 

-4E( 4J ) er ,xr, ^(f^n^'M^ 

- 2 D(;j)er 1 xr 2 1 [v*^<0] ll- 2 ^ ~ z ij\\ 2 + CS( l , :) )er 1 xr 2 (%) 2 ' 
where, C is a constant, yT\ = max(— yy, 0) and — max(j/y,0). This inequality follows easily from the fact 
that, for any G T 1 x T 2 , 

(i) f OT (s,A^,(2; fei )( M)e rixr 2 ,%-) < /^" +1 ^(s, X* a >*, (y k i) {k , l)er i xr 2 , *«), 

(ii) For any (ttfcOf^Oer 1 xr 2 G such that My = and Uki > 0, for (k,l) ^ 

f ^ m (s, Xt' x , (y kl + u kl ) (iJ)er i x r*, ztj) > f ij ^ m (s, X* a >*, (y k i) (k , l)e ^ x ^,Zij)- 

(iii) fij depends only on Zij and not on the other components z k i, [k, I) ^ (i, j). 
Consequently, we have 

P - a.s., Y ij > n ' m < Y ij ' n+1 ' m , G r 1 x r 2 . 

In the same way one can show that Y ,3 ' n ' m+ < Y l3 ' n,m . Finally, the inequalities of (|3.4p are obtained by taking 
s = t in d3HD in view of the representation (UTS) of Y i ^ n ' m by v ij > n < m and X l ' x . U 

We now suggest two approximation schemes obtained from the sequence (Y l3 > n > m , (i,j) G T 1 x r 2 )„ jm of the 
solution of the system (13.11) . The first scheme is a sequence of decreasing reflected BSDEs with interconnected 
lower obstacles and the second one is an increasing sequence of reflected BSDEs with interconnected upper 
obstacles. 

Let us first introduce the decreasing approximation scheme by considering the following system of reflected 
BSDEs with interconnected obstacles: V(i,j) G T 1 x T 2 , 

yij.m g S 2,l^ zij,m g y2,d and ^ij,m g j2,+ . 

yy,m = h ij ( X t,x^ + jT pj,m^ x t,x f (Y r kl ' m ) (k iJ)er i x r 2 , Zp' m )dr + Rf m - Kf m - jf Z^ m dB r , 
' Y^ m > max {Y^ m - gjs, X**)}, s < T, (3 ' 5) 

k ff(Yp m - nua i6(P) -»{yW- m - gjs, Xr;\ulls:>"> = 0, 
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where, for any G T 1 x T 2 , m > and (s,y,z i: >), 

fV' m (s,Xi'*,y, Z v) :=/ ij '(*,^' x ,(y* , )(*,i)en x r»,2 y )-m(y«- min + g^s, X^))) + . 

ie(r 2 )-J 

Thanks to Assumptions (H1)-(H3) and (|2.3|l . by Theorem 3.5 in Hamadene and Zhang (2010) ([17]), the solution 
of (13.51) exists and is unique. In fact, this holds again under weaker assumptions (see Hamadene and Morlais 
(2011)f[T9"]. Theorem 1). For sake of completeness, a statement of this recent result is given in Appendix (A2). 
Moreover, we have the following properties. 

Proposition 3.2. For any £ T 1 x T 2 and m > 0, we have 
ft) 

E[ sup \Y* j ' n ' m -YV> m \ 2 } -> 0, as n -> oo, (3.6) 

t<s<T 

fti) 

P-a.s., Y ij ' m >Y ij ' m+x . 

(Hi) There exists a unique A-uplet of deterministic continuous functions (v' m ) fatter 1 xr 2 * n n g such that, for 
every t < T, 

yij,m = x**), s e [t, T]. (3.7) 

Moreover, for any G T 1 x T 2 and (t,x) G [0,T] x R k , v ll - m (t,x) > v i i< m + 1 (t, x) . 

Finally, (^'"Jiijjepxr 2 is the unique viscosity solution in the class Tl g of the following system of variational 
inequalities with inter- connected obstacles. V G T 1 x T 2 , 

mm{v lj ' m (t, x) — max (#-?' m (i, x) — g (t, x)); 

feeCr 1 )-^ — lk 

-d t v^ m (t,x) - b(t,x)D x v ij ' m (t,x) - ±Tr{o-o- T (t,x)Dl x v i i> m (t,x)) ^ g ^ 

-f^ m (t, x, (v kl > m (t, a;)) (fc)06 ri , <? T (t, x)D x vV> m (t, x))} = 0, 

v ij - m {T,x) = h l i{x) 

where, 

f J ' m (s,x, (y fe ')(fe,/)erixr2,2;) = f v {s,x, {y kl )fa l)er i xr 2, z) - m(y v - min (y d + g jt (s, x))) + . 

/e(r 2 )-3 

Proof. Let us prove (i). For this, it is enough to consider the case m = 0, and we will do so, since for any 
G T 1 x T 2 , the function 

(s,x,(y kl ) {u)er i xr 2) — > -m(y ij - min (y a + g 5l {s, x))) + 

;e(r 2 )-J 

has the same properties as f 1 ^ displayed in (H1)-(H2). 
To begin with, let us show that for any i,j and n > 0, 

P-a.s., Y ij ' n '° < Y ij '°. (3.9) 

First andt w.l.o.g. we may assume that / IJ is non-decreasing w.r.t. (y kl )(k,i)er 1 xr 2 , since thanks to Assumption 
(H2), it is enough to multiply the solutions by e ro ', where w is appropriately chosen in order to fall in this latter 
case, since / y is Lipschitz in y l K Now, for fixed n, let us define recursively the sequence ' ,J '' n )fc>o as follows: 
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For k = and any G T 1 x T 2 , we set F°>«'« : = Y i3 >° and, for any fc > 1, let us define (Y k ' ij ' n , Z k > i3 < n ) G 
5 2,i x y2,d ag the so i ut i on of the following system of BSDEs: V(i, j) G T 1 x T 2 , 

_ d y*,ii,n = fV^x**, (Yt 1 ' pq ' n )( P , q )erixr^Z^> n )ds+ 

n{Y s k >»> n - ma^i^yiiY*- 1 ' 1 ^ - g. k (s,Xl> x ))}-ds - Z^> n dB s , (3.10) 

The solution of (I3.10[) exists since it is a multi-dimensional standard BSDE with a Lipschitz coefficient, noting 
that (Y k > pq ' n )( p ^er 1 xr 2 is already given. Since, n is fixed and the coefficient 

^(^(/Voer^r*,^ max ,(y«- s (a,^)}" 

/e(r 1 )- 1 _ " 

is Lipschitz w.r.t. ((y kl )(k,i)er 1 xr 2 ^ z^)> the sequence (Y k,l3 ' n )k converges in S 2,1 to y^' n >°, as k — > oo, for any 
i,j and n. 

Using an induction argument w.r.t. k, we prove that, for any i,j and n, 

P - a.s., Y k ' ij ' n < Y fc > 0. 

Indeed, for fc = the inequalities hold true through the definition of Y 0,l3,n . Assume now that these inequalities 
are valid for some k — 1, i.e., for any i, j and n, 

P - a.s., y*-i.w» < y «.o. (3.11) 



Thus, taking into account that Y 13 ' satisfies (|3.5[) with m = and the fact that f 13 is non-decreasing w.r.t. 



(y fe 0(fe,()er 1 xr 2 then for any n, it holds 



<fH S ,X^,(Yp-°) (p . q)erlxr 2,z^). 

Using the standard comparison result of solutions of one dimensional BSDEs we obtain that, for any n,i,j, 
Y k ' l i' n < F'i" . a.s. Thus the property (|3.11[) is valid for any k. Taking the limit as k tends to oo, we obtain 

Jan. 



We can now apply Peng's monotonic limit (see Peng (1999) ([H])) to the increasing sequence (Y y '"' ) n >o- This 
yields the existence of processes Y ij , Z %3 and K ij , G T 1 x T 2 , such that: 

(a) Y %3 is RCLL and uniformly P-square integrable. Moreover, for any stopping time r, lim n ^oo /*■ Y* 3 > n '° = Y* 3 . 

(b) K %3 is RCLL non-decreasing, Kq = and for any stopping time r, lim„_ J . 00 K l T 3 ' n '° = K %3 , P — a.s. 

(c) Z ij G n 2 ' d and for any p G [1, 2), 

limE[/ \Zi j ' nfi -Zi j \ p ds] = 0. 



(d) For any G T 1 x T 2 , the triple (Y ij , Z ij , K ij ) satisfies: Vs < T, 



Y? = h i3 (X^) + jf /«(r, X*>*, (Y r k % jl)er r xr2 ,ZV)dr + K% - K» - £ Z»dB r 
{ Y? > max fc6(rl) - 4 {^ -g. k (s,Xl>*)}. 

The remaining of the proof is the same as in Hamadene and Zhang (2010) (([U]), Theorem 3.2) and it mainly 
consists in proving both the continuity of Y %3 and K %3 and the Skorohod condition (these properties are deduced 
using the no- free loop property (H4)). Finally we obtain that the triple (Y %3 , Z %3 , K 13 ) satisfies (|3.5[) with m = 0. 
Hence, by uniqueness of the solutions of (|3.5I) . Y 13 = F 1 ^ a.s, which completes the proof of (i). 
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(ii) is an immediate consequence of (i) and Proposition 13. II 



We now establish (iii). The existence of the deterministic continuous functions (w y ' m )(i,j)er 1 xr 2 such that for 
any s £ [t,T], Y^' m = v l ^ m {s,X s ), £ T 1 x V 2 , m > 0, and which satisfy the system fl-£HJ|, is obtained in 

[19] (see Appendix A2, Theorem^. Finally as Y ij ' m > Y^' m+1 a .s., we deduce that v l ^ m > v ij ' m+1 , taking 
into account (ii), which completes the proof. M 

We now introduce the increasing approximation scheme by considering the following system of reflected BSDEs 
with interconnected obstacles: for any (i,j) £ T 1 x T 2 , 

' Y ij ' n £ S 2 , Z ij ' n £ H 2 ' d and K ij ' n £ A 2 '+ ; 

rr = h^(x^) + j;'r'i, a;-. <x?> n ) w)& i xF *,z?' n )dr - m> n -Kf n ) - £z?rdB r , 

rJ^K^mmjr^+g^Xln}, s<T, (3.13) 

T 



where, for any (i,j) £ T 1 x T 2 , n > and (s, y, z % i 



r' n (s,X^,y,zV) = ri(s,xr,(y kl ) (k , l)erlxr2 ,zy) + n(yv ~ max .{y« - g (a, X^)}) . 

Thanks to assumptions (H1)-(H3) and (H4)-([%3|). by Theorems 3.2 and 3.5 in Hamadene and Zhang (2010) (fi7\) 
the solution of p.5[) exists and is unique. 

Below is the analogous of Proposition ^. 21 We do not give its proof since it is deduced from this latter proposition 
in considering the equation satisfied by (— Y? 3,n , —Z l3 ' n ,K VJ,n ). 

Proposition 3.3. For any £ T 1 x T 2 we have 
(<) 

E[ sup \ Y l ^ n ' m -Yf n \ 2 ] ->0 as m -> oo. (3.14) 

t<s<T 

(«i) For an?/ n > 0, 

p- a .s., r''" • !_"•'••'. 

(iii) For every n > 0, there exists a unique K-uplet of deterministic continuous functions (w fei '")(fc,z)er 1 xr 2 ™ n g 
such that, for every t < T , 

Yf n = V ij ' n {s,Xl' x ), t<s<T. (3.15) 
Moreover, for any n > and £ T 1 x T 2 , 

3i y>n (i,aj) < v ij ' n+1 (t,x), (t,x) £ [0,T] x M fc . 

Finally, (u" J ' n Wj)er 1 xr 2 *s i/ie unique viscosity solution in the class H g of the following system of variational 
inequalities with inter- connected obstacles: V(i,j) £ T 1 x T 2 , 

max-fv'^fi, x) — min (v d > n (t, x) + gi k (t, x)); 
ie(r 2 )-j 

-d t v^ n (t,x) - b{t,x)D x v^ n (t,x) - ±Tr(aa T {t,x)D 2 xx v^ n (t,x)) (3 lg) 

-/ lJtn (t, x, (v kl > n (t, x)) {k , l)er i xr 2,a T (t, x)D x v^ n (t, x))} = 0, 

k v ij ' n (T,x) = h i3 (x) 



where, 



f ij ' n (s,x,(y kl ) (k>l)er x xr 2,z) = f l 3 {s,x,{y kl ) {ktl)er i xr 2,z) + n{y ij - max (y kj - g (s, x))}~ 
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For (t,x) G [0,T] xl'e [0,T] x M' c and (ij) G T 1 x T 2 , let us define 

#(i,x):= lim ^'' m (i,a;) and v ij (t,x):= \im v ij ' n (t,x). 
Then, as a by-product of Propositions 13.21 and 13.31 we have 



Corollary 1. For any G T 1 x T 2 , the function v 13 (resp. ) is use (resp. Isc). Moreover, v 1 ^ and v 13 

belong to H g i.e., there exist deterministic constants C > and 7 > such that, for any (t,x) G [0, T] x M. k , 

\v ij (t,x)\ + \v l3 {t,x)\ < C(l + |a;| 7 ), (t,x) G [0,T] x M fe . 
Finally, for any (t,x) G [0,T] x M fe . 

^'(£,2) < u ij '(t,x). 

Proof. For any (i,j) G T 1 x T 2 , v 1 ^ (resp. v^) is obtained as a decreasing (resp. increasing) limit of 
continuous functions. Therefore, it is use (resp. Isc). Next, for any and n,m, 

u «J>,m(t J a;) < v ij ' nfi (t,x), (t,x) G [0,T] x M fe , 

as the sequence (v %J ' n ' m ) m >o is decreasing. Thus, taking the limit as m —> 00 we obtain 

Now, using (HU and dH) it follows that, for any i < T and s G [t,T], Y* j < n <° = v ij ' n '°(s,Xj' x ) and the processes 
yy',n,o conver g es in 5 2 , as n — ► 00, to Y l ^° solution of (|3.5p with m = 0. Furthermore, by (|3. 71) . there exists a 
deterministic continuous function u 4 - 7,0 of class II g such that for any t < T and s G [t, T] , F^' = v 1; >< (s, Xl' x ). 
Then, taking s — t and the limit as n — > 00 to obtain 

v ij (t,x) := lim v ij ' n (t,x) < lim ac) = « y '°(t,a:), V(*,£c) G [0,T] x R fc . 

But, -D^' and v ij ' n belong to II g and v ij ' n < y_ ij > n+1 . Thus, %f° belongs to for any G T 1 x T 2 . In the 
same way one can show that v 1 ^ belongs to U g , for any (i,j) G T 1 x T 2 . The last inequality follows from p.4j) 
and the definitions of w y and . M 

4 A comparison result 

In this section we investigate some qualitative properties of viscosity solutions of the system (12.11) . In particular, 
we show in Corollary [5] below that if the system (|2.1j) admits a viscosity solution in the class II g , then it is unique 
and continuous. 

We first show that if (u y )( ij ) g r 1 xr 2 (resp. (ifl^Vjjjer^xr 2 ) is & n u $ c subsolution (resp. Isc supersolution) of 
(|2.ip which belongs to II g , then for any (i,j) G T 1 x T 2 , < . To begin with, we give an intermediary result 
which is required in the second step of the proof of the comparison result. 

Lemma 4.1. Let u := (u i:l )u t j\^r 1 xT 2 (resp. w :— (w u )(i. :) )er 1 xr 2 ) oe an usc subsolution (resp. Isc supersolu- 
tion) of \2.1\) . Let (t,x) G [0,T] x R k and let f(t,x) be the following set. 

f(t,x) := G T 1 x r 2 ,u ij (t,x)-w i ' j (t,x) = max ^{u kl (t,x) - w kl (t,x)}. 

Then, there exists («o>io) € T(t,x) such that 

u iojo (t, x) > L iojo [u\(t,x) and w iojo (i, x) < U iol ° [w](t, x). (4.1) 
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Proof: Let (t, x) G [0, T] x R fe be fixed. Since the set L 1 x L 2 is finite then f (t, x) is not empty. Next, let us show 
the claim by contradiction. So for any G T(t,x) either u lJ (t,x) < L 13 [u](t, x) or w l3 (t,x) > U t3 [w](t,x) 

holds. Let us first assume that: 

u' 3 (t,x) < L ij [u](t,x). (4.2) 
Then, there exists some k G (L 1 ) \ such that 

u ij (t,x) < L»[u\(t,x) = u k i(t,x)-g. k (t,x). 

But, since w is a supersolution to (|2.1[) we also deduce 

which implies that 

u 13 (t, x) - u kj (t, x) < -g. k (t, x) < w ij (t, x) - w kj (t, x) 

and then 

u ij (t,x)- w ij (t,x)< u k3 {t,x)~ w kj (t,x). 

Since, by assumption 6 T(t,x), the two previous inequalities are instead equalities, (k, j) belongs to T(t,x), 
k 7^ i and finally it holds that 

u« (t, as) - (t , x) = -g ik (t,x) = w» (t,x) ~ w" 3 (t,x). (4.3) 

Next, if (|4.2[) does not hold, then necessarily u lJ (t,x) > L lJ [u](t,x) and w %3 (t,x) > U l3 [w](t,x). As (u l3 )i.j is a 
subsolution of (|2.1| we obtain 

u ij {t,x) < U rj [u]{t,x) < u a (t,x) + gji(t,x), I 6 L 2 " J '. 

On the other hand, for some index I G 3 it holds 

(i, x) - w il (t, x) > g-ji {t,x)> u 13 (t , x) ~ u a (t , x) . 

Once more as (i,j) G T(t,x) then the previous inequalities yield that (i,l) G T(t, x), I ^ j and 

u ij {t,x) ~u l \t,x) =gji(t,x) =w l3 (t,x) -w il (t,x). (4.4) 

Repeating now this reasoning as many times as necessary, and since T 1 x T 2 is finite, there exits a loop 
(ii, ji), . . . , (in-i, Jn-i), {tN, 3n) = («i, ji) such that 

g=l,jV-l 

which contradicts Assumption (H4), whence the claim is proved. M 
We now give the main result of this subsection. 

Theorem 4.2. Assume that u = (u 13 )(i.j)er 1 xr 2 (resp. w = (w %3 )uj)^t 1 xt 2 ) is an use (resp. Isc) subsolution 
(resp. supersolution) of the system \2.1\) such that, for any G L 1 x L 2 , both u 13 and w 13 belong to Yl g , i.e., 

there exist two constants 7 and C such that 

G L 1 x L 2 , V(t,x) G [0,T] x R fc , \u l3 (t,x)\ + \w ij (t,x)\ < (7(1 + M 7 ). (4.5) 

Then, it holds that for any (i,j) G L 1 x T 2 , 

u lj {t,x) < w ij (t,x), (t,x) G [0,T] x M fc . 
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Proof. Let us proceed by contradiction and let (u lJ )(i,j)er 1 xr 2 (resp. w — (w l:l )ii i j\ e r 1 xT 2 ) be use (resp. Isc) 
and a subsolution (resp. a supersolution) of the system (|2.1[) such that there exists £o > and (to, xo) € [0, T] xM. k 
such that 

:((u ij -w ij )(t ,x ))>e . (4.6) 



max i 



Next, w.l.o.g. we may assume that for any G F 1 x F 2 , 

lim {u ij - w ij ){t,x) = -co. (4.7) 

\x\ — >og 

Indeed, if this is not the case, one may replace w ij with w 13 ' defined by 

w ij ^^(t,x) =w ij (t,x)+'&e-< lt \x\ 2 -i +2 ,(t,x) G [0,T] x R k , 

which is still an use supersolution of (|2.1[) for i9 > and fx > hq which satisfies (|4.7|l (a proof of the supersolution 
property for good choices of <d and fi can be found in e.g. Pham (2009) ((|25. ), pp.76). Therefore, it suffices 
to show that u t3 (t,x) < w 13 ' '^(t, x), (t,x) G [0,T] x K fe , since, by taking the limit as — > 0, we deduce that 
U ij (t,x) < w ij (t,x), (t, as) G [0,T] xl fe . Thus, assume that gl)]) and (@2J are satisfied. Then, there exists R > 
such that 

max (M)e[0 , T]xKfe maxy^u' 3 - u> l: 0(i, x)} = max (M)E[0iT | x , J ( j ) max;j{(u lJ - x)} 

= majc,^!!^ - w ij ')(t*,x*) > eo > 0, 

where, (**,£*) G [0,T) x B(0,i?), where, B(0,R) denotes the ball in with center the origin and radius R, 
since by definition u ij (T,x) < w i3 (T,x), for all G T 1 x T 2 . 

The remaining of the proof is obtained in two steps: the first step which is the main one establishes the comparison 
result under the additional condition (|4.9[) and the second step provides the proof in the general case. 



Step 1. Let us make the following assumption on the functions (/ y ')(t jjer^xr 2 - F° r a U 

A > c(/' iJ )(A - 1), (i,j) er 1 xr 2 , (t,x,y,z) G [0,T] x R k x R A+d , and (u, u) G M 2 s.t. u>v, 

rHt,x,[y-^\u], z) - f 3 (t,x,[y- {t3 \v], z) < -X(u-v), (4.9) 

where, c(f ij ) is the Lipschitz constant of f lJ w.r.t. {y kl )(k.i)er 1 xr 2 • Next, let (io,jo) be an element of T(t*,x*) 
that satisfies (I4.1[) . For n > 1, let $^ 0j0 be the function defined as follows. 



(t,a;,!/) := (ti^^a;)-^"^!/))-^.^!,!/), G [0,T] x 



where, 



y) := n\x - y\^+ 2 + \x- x*\ 2 ^+ 2 + (t - t*) 2 . 
Since ^ i ° ,3 ° is use in (t,x,y), there exists (t n ,x n ,y n ) G [0, T] x B(0,R) 2 such that 



Moreover 



J A (t* , x* , x*) = ■*> (t* , ) - u> >*> (f ,i*)< J" (i* , x*) - «;*> >*> (t*,x*) + K {t n ,x n ,y n ) 
<u l °^(t n ,x n )-w l °^{t n ,y n ). 



13 



The definition of 4> n together with the growth condition of u u and it/ 1 -? implies that (x n — y n ) n >i converges to 0. 
Next, for any subsequence {{t ni ,x ni ,y ni ))i which converges to (i, x, x) we deduce from (|4.10p that 

u iojo (t*,x*) - w iojo (t*,x*) < u io ' jo (t,x) - w iojo (t,x), 

since u l °3" is use and w io ^° is Isc. As the maximum of u SoJO — w 10 ^ on [0, T] x B(0, R) is reached in (t* ,x*), then 
this last inequality is actually an equality. Using the definition of <p n and (|4.10p . we deduce that the sequence 
((t nj x ni y n ))n converges to (t*,x*,x*) which also implies 

n\x n -y n \ 2l+2 -> and (u iojo (t n ,x n ),w iojo (t n ,y n )) -> (u ioj ° (f, x*), w ioj ° (t*, x*)), 

as n —> oo. This latter convergence holds, since we first obtain from (14. 10)) and in taking into account that u Zo: > 
and w l ° 30 are Isc and use respectively, 

u ^(t*,x*) - w l "^{t*,x*) < ljm n (u^>{t n ,x n ) - w l ^{t n ,y n )) < Ihn"„(u^° (t n , x n ) - w^"{t ni y n j) 

<mi n u iojo (t n ,x n )-]xmnW iojo (t n) y n )) < u l °3°(t*,x*)-w^°(t*,x*). 

Thus the sequence (u l " J "(t n ,x n ) - w lojo (t n , y n )) n >o converges to u loJO (t* , x*) - w loJO (t* ,x*) as n — > oo. Then 

]hn n u lo: > (t n ,x n ) = u iojo (t* ,x*)- w la 3° (t*,x*) + lim„ w iojo (t„ . y n ) > u iojo (t* , x*) > lmT„ u loJ0 i„, x„). 

It follows that the sequence {u 10 ^ (t ni x n )) n converges to li 10 -* (t* , x*) and then (tt! loJO (t n ,y n )) n converges also to 
w io3o {t*,x*). 

Next, recalling that u % °^° (resp. w % °i°) is use (resp. isc) and satisfies (|4.ip . then, for n large enough and at least 
for a subsequence which we still index by n, we obtain 

u lo3o {t n ,x n ) > max (u kjo (t n ,x n ) - g iok (t n ,x n )), (4.11) 
fcG(ri)- 

and 

w ioj °(t n ,x n ) < min (w iol (t n ,x n ) - g jo i(t n ,x n )). (4.12) 
/e(r 2 )-^ 

Applying now Crandall-Ishii-Lions's Lemma (see e.g. [5] or [10], pp.216) with <&j£ j0 and at the point (t n , x n , y n ) 
(for n large enough in such a way that this latter triple will belong to [0, T] x B(0, R) 2 ), there exist (p™, g™, M") € 
J 2 ' + (u ioJo )(t n ,a; n ) and (p£,g",M«) e J 2 '" (w* ' )(*»,!/.,) such that 

K - PS = 9 t ^ n {t n ,x n) y n ) = 2(t„ - t*), g™ (resp. g") = d x ip n {t n ,x n ,y n ) (resp. - d y ip n (t n ,x n ,y„)) and 

f M " ° Wn + ^, (4.13) 

where, A n = D 2 x y ^ip n (t n , x n ,y n ). Taking into account that {u l ^)i.j (resp. (w* J )t,j) is a subsolution (resp. 
supersolution) of (I2.ip and using once more (|4.11l) and (|4.12l) we get 

- V l - b(t n , x n ) T .ql - ^Tr[(aa T )(t n , x n )M^\ - (t n , x n , (u ij (t n , x n )) (lj)er i xr 2,a(t n ,x n ) T .q^) < 0, 

and 

-pi - b{t n ,y n ) T .ql - ^Tr[(aa T )(t n ,y n )M^] - f ojo (t n ,y n , (w ij (t n ,y n )) {itj)er i xT 2,a(t n ,y n ) T .q%) > 0. 

Taking the difference between these two inequalities yields 

-(Pu-Pl) - (b(t n ,x n ) T .q2 - b(t n ,y n ) T .C) - ±Tr[{aa T (t n ,x n )MZ - aa T (t n ,y n )M%}} 

-{r o3o (tn,x n ,(u^(t n ,x n )) (Mj)er i xr 2,a(t n ,x n ) T • (4.14) 
-P 03a (tn,ynAw t Ht n ,y n )) {lJ)er i xr 2,a(t n ,y n ) T .q™)} < 0, 



14 



and then 

-{P' yJa (t n , x n , {u 13 (t n , a; n ))(i,j)eri , <r(i n , a;„) T .<7™) 

-f %ojo (tn,y n , (w y (^,yn))( l j)erixr2,cr(i Il ,y„) T .q^,)} < £„, 

with limn^oo p„ < 0, using the fact that all the terms in the first line of (I4.14p are converging sequences. 
Linearizing f 1 " 30 (see Appendix A3), which is Lipschitz w.r.t. (y IJ )(jj) er i xr 2, and using Assumption (|4.9[) , we 
obtain 

A^'"^,^) -w iaja (t n ,y n )) - ^ 0«(u y (t n ,x n ) - w ij (t n ,y n )) < g n> 

(ij)enxr 2 ,(i,i)^(i 3o) 

where, is the increment rate of p 030 w.r.t. y 13 which is uniformly bounded (w.r.t. n) and non-negative 
thanks to the monotonicity assumption (H2). Therefore, 

\{u l030 {t n ,x n ) -w io3 °(t n ,y n )) 

< E(i,i)grixr»,«j)#« jo) ®ii j {u ij (tn,x n ) - w i3 {t ni y n )) + g n 

< c(f iojo ) x T,(ij)er^xr^(i,j)^i ,j )(( uij ( t n,x n ) - w l3 (t n ,y n ))+ + g n . 

Taking the limit as n — >• oo we obtain 

X(u iojo (t*,x*) -w iojo (t*,y*)) 

< IW^oo c(/ iM '°) (E(i,j)erixr=,(i,j)^(io,jo)( uiJ (*«' a: «) _ "^'C^l/n))" 1 ") 

< c(/ W0 ) fe (tj)erlxr2 (jJ . )#(lo Jo) IW rl ^ 00 (^(i„,x„) - w <J '(i„,j/„)) + ) 

since u 4 - 7 (resp. it; 4 - 7 ) is use (resp. isc). As (io,jo) belongs to r(i*,a;*), we obtain 

A(u Ioio (r,a;*) - w 4030 (t*,y*)) < c(/ Wo ) ((A - l)(u iojo (t*,x*) - w io3 ° (t* , x*))) . 
But this is contradictory with and (|4.9p . Thus, for any G T 1 x T 2 , it 4 - 7 < w 4 - 7 . □ 

Step 2. The general case. 

For arbitrary A € R, if (u^Wjjer 1 xr^er 1 xr 2 (resp. (w y )(i,j)er 1 xr 2 ) is a subsolution (resp. supersolution) of 
(|2.ip then u* 3 (t, x) = e xt u 13 (t, x) and w ij (t, x) = e xt w 13 [t, x) is a subsolution (resp. supersolution) of the following 
system of variational inequalities with oblique reflection. For every G T 1 x T 2 , 

minj-D 4 - 7 (t, x) - max le ^-i{v l3 (t,x) - e xt g u (t 7 x)}; 



ma&{v %3 {t,x) - rRm ke(r2 )-] {e xt g jk {t, x) + v 13 (t,x)}; -d t v 13 (t, x) ~ Cv 13 (t, x) + Xv zj (t,x) 
-e xt f 3 (t, x, {e~ xt V 3 {t, x)) (iJ)&r i xr2 , e~ xt a T (t, x).D x v»(t, x))} = 0, 
v ij (T, x) = e XT h i3 (x). 

But, by choosing A large enough the functions 

F i3 (t,x,(u k % il)er ^^,z) = -Xu i3 +e xt r 3 (t,x,(e- xt u kl ) {k ^ rlxT2l e- xt z), G T 1 x T 2 , 

satisfy Condition (|4.9j) . Hence, thanks to the result stated in Step 1, we have u 13 < v %3 , G T 1 x T 2 . Thus, 
u %3 < v K for any G T 1 x T 2 , which is the desired result. ■ 

As a consequence of this comparison result, we obtain the following one related to uniqueness of the solution of 

USED- 
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Corollary 2. If the system i2.1\) admits a viscosity solution in the class Ii g , then it is unique and continuous. 

Proof. Indeed, assume that ^^u^^r^-xT 3 1S a solution of (|2.1j) that belongs to II g . Then, thanks to 
the previous comparison result, for any G T 1 x V 2 we have v 1 ^* < vl 3 . Thus, v t3 '* — vl 3 and then 

v %3 — vl 3 — v %3 '*, which means that v 13 is continuous. Next, if {v l3 )u,j)&r 1 xv 3 is another solution of (|2.1|) in the 
class n s , then it is also continuous and by the Comparison Theorem 14.21 we have v 13 < v 13 and o' 3 > v %3 , Hence, 
v *j — € r 1 x r 2 , i.e., uniqueness of the solution of (|2.1I) . ■ 



5 Viscosity solution of the system (12.11) 



In this section we prove that the family {v lJ )i^ constructed in Section 3 provides the unique continuous solution 
in viscosity sense of the system (|2.1[) . For sake of clarity, the proof is divided into several steps. 

Proposition 5.1. The family {v tJ )(i j)er 1 xr 2 is a viscosity subsolution of the system \2.1\) . 

Proof. First recall that for each G T 1 x T 2 , v %3 is use, since v %3 = lim m \ v l3 ' m , where v lJ ' m is a 

continuous function solution of the system (|3.8[) . Thus, for any (i,j) G T 1 x T 2 , it holds that v tJ '* = v %0 . Next, 
at T we have 

v 13 (T, x) — lim \ v ij ' m (T, x) = h ij (x) , i6l'. 

m 

We shall now prove that, for any (t,x) in [0,T) x IR fc , any G T 1 x T 2 and (p,q,M_) in J + v tJ '(t, x), 



(5.1) 



min[(^ - L^')(i,x),max{(^ - U i3 ){t,x), 

-p-b(t,x)q- \Tv{aa T {t,x)M) - f 3 ' (t,x,{v kl (t,x)) (kil)er i x r*,o- T (t,x).q)}\ < 0, 

with Z 1 - 7 and L™ defined as follows: 

L i3 (t,x) = ^ max (v k3 (t,x) - g. k {t,x)) and U ij (t,x) = ^ mm (v il (t,x) + gji(t,x)). 

Now, let G r 1 x r 2 be fixed. Then it is equivalent to show that, either 

(B* 3 "-Z y )(t,x)<0, (5.2) 

or 

max{(# - U ij ){t,x),-p-b{t,x)q-^TT(aa T {t,x)M) - P '{t, x, v(t, x), a T {t, x).q)} < 0. (5.3) 

If (|5.2p is satisfied then the subsolution property (|5.1|) holds. Therefore, from now on, we suppose that there 
exists eo > such that 

v ij {t,x) >L ij (t,x) + e , (5.4) 

and show (|5.3[) . Thanks to the decreasing convergence of (u y ' TO ) TO >o to v i3 , G T 1 x T 2 , there exists mo such 
that for any m > mo, we have 

v ij ' m (t,x) > L ij ' m {t,x) + (5.5) 
Next, by continuity of v 13 ' m and I/ 3 ' m , there exists a neighborhood O m of (t,x) such that 

^ m (i', x 1 ) > L^ m (t', x 1 ) + ^, (*', x') G m . (5.6) 

Now, by Lemma 6.1 in [5] there exists a subsequence ((tk,Xk))k>o such that 

(tfc.aifc) -^k^oo (t,x) and ^(£,2;) = lim v ij ' k (t k ,x k ). (5.7) 
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Moreover, there exists a sequence (p k ,q k ,M k ) G J + (u t5 '' k (£fe, such that 

lim (pk,qk,M k ) = (p,q,M)- (5-8) 

re— foo 

* > i ' . the subsequence ((ife, ^fe))A:>o can be chosen in such a way that for any k > 0, (t k ,x k ) G 0/c- Applying now 
the viscosity sub-solution property of v^' k (which satisfies (|3.8p ) at (t k ,x k ) and taking into account of (|5.6[) we 
obtain 

~Pk-b{t k ,x k ) T .q k - ^Tr(aa T (t k ,x k )M k ) - f j ' k (t k ,x k ,(v pq ' k (t k ,x k )) {Ptq)er ^xr^cr(tk,x k ) T q k ) < 0, (5.9) 
where, once more, 

f 3 ' k (s,x, {y pq ) [pA) ^xT^z) = f j (s 7 x } {y pq ) Mer i xr 2,z) - k(y ij - min [y a + g 3 i{s, x))) + . 

ie(r 2 )-3 

Next, thanks to the boundedness of the sequence ((t k ,x k )) k >o, the uniform polynomial growth of v pq ' k k > 0, 
(by Proposition 13.21 and Corollary [lj, the assumptions (H0)-(H2) on b, a and /*•?', and the convergence of 
((p k ,qk, M k )) k to (p, q, M), we deduce from (|5.9II that 

e k := (v l3 ' k (t k ,x k ) - min (v U ' k (t k ,x k ) + g~ji(t k ,x k ))) + -> 0, fc -> oo. 
(s(r 2 ) j 

But, for any fixed (t,x) and fco, the sequence (v ll ' k (t,x)) k > ko is decreasing and then for any k > ko > too, 
v ij ' k (t k ,x k ) < min (v il ' k (t k ,x k ) +g jl (t k ,x k )) + e k < min (v il ' ko {t k , x k ) + g 3 i(t kl x k )) + e k . 

Taking now the limit as k +oo, in view of the continuity of v ' °, we get 

limv^ k (t k ,x k ) = v ij (t,x) < mm(v il ' ko (t,x)+g :jl (t,x)). 

k j^l 

Finally, passing to the limit as fco goes to +oo to obtain 

v ij {t,x) < mm(v a {t,x)+g jl (t,x)) = U l °{t,x). (5.10) 

Let us now consider a subsequence of (k), which we denote by (ki), such that for any (p, q) G T 1 x T 2 , the sequence 
(v pq ' k (t kl ,x kl ))i is convergent. This subsequence exists since the functions v % ^ kl are uniformly of polynomial 
growth (by Proposition 13.21 and Corollary [T]) . Then, taking the limit w.r.t. I in equation (|5.9[) , we obtain 

-p-qb(t,x) - |Tr(crcr T (i,x)M) < lim;-^ /' y < fei (t kl , x k , (v pqM (t kl , x kl ))( p , q)eV i xT 2 , a(t kl , x kl ) T q kl ) 

< lim/^oc f ZJ (t kl) x k , (v pqM (t kl ,x kl )) {p ^ q)er i xr 2,<j(t kl ,x kl ) T q kl ) (5.11) 
= P J {t, x, (lim^oo v pq - kl (t kl , x kl ))(p, g ) e r 1 xr 2 , ^) T <z), 

since / lJ is continuous in (i, £, y, z). Now for any (p, q) G T 1 x T 2 , since -D P9 '™ is continuous and v pq,n > i; pg '™ +1 , 
Vn > 0, it holds that 

v pq '*[t,x) = v pq {t,x) = lim" v pq ' n {t',x'), (t,x) G [0,T] x R*. 

Therefore, for any (p, g) G T 1 x T 2 ((p, g) ^ we have 

«M(t,aO > lim v pq ' k '(t kl ,x kl ) and = lim v ij > k ' (t kl ,x kl ). (5.12) 
As / y is non-decreasing w.r.t. ?/ fci , (k,l) G T 1 x T 2 , (k,l) ^ we deduce from (15. lip and (|5.12[) that 

-p-qb(t,x) - ±Tr(aa T (t,x)M) < f i(t, z, (v pq (t, x)) Mer i xr ,, a(t, x) T q). (5.13) 
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Finally, under the condition (|5.4I) . the relations (|5.13p . (|5 . 10[) imply that f|5 . 3[) is satisfied. Thus v l i is a viscosity 
subsolution for 



- L ij ) (t, i),max{ (v ij - U l3 )(t, x) 

r £- \T±{oo T {t,x)M) - i'- I ' • ! i i, 

v ij (T,x) = h ij (x) 



p - b(t, x) T q - ±Tr(aa T (t, x)M) - /« (t, x, (v w (t, x)) (fe , )erl xr , , a T (t, x).q)}} = 0, 



Since is arbitrary, jg^xr 3 is a viscosity subsolution for (|2.1|) . This finishes the proof. ■ 

Proposition 5.2. Let too fee /ixed in N. Then, the family (v lJ ' mo )fi t j-)^r 1 xT 2 ^ a viscosity super solution of the 
system \2.1\l . 

Proof. First and thanks to Proposition O we know that the triples ((Y^' m °, Z^' m °, K l ^ ma )) {l ^ )eT i xr2 
introduced in Q3.5p . solve the following system of reflected BSDEs: For every 6 T 1 x T 2 , 

y ij.mo g 52,^ ^y,m g %2,ei and ^y,m g ^2,+ . 

Y i3,m =hv(X% x ) + jj f l 3> mo (r, X^ x , (Y r kl < mo ) {kJ)er i xT 2 , Z^ mo )dr + Rf m ° - Kf m * - Z*?' m °dB r 
Y?™ > ^JY*™ -g. k (s,X^)}, s<T; 

{ Jo(Y s ij ' m ° -maxfcg^)-^?^ - g. k {s,X^)})dRf^ = 

(5.14) 

where, for any 6 T 1 x T 2 and (s,y, z 4 - 7 ), 

f^( S ,X^,y,^) = /* J >,X^(y%, 0er i x r>,** J ') -m (^ - min .(/ 

;g(r 2 )- 3 

Furthermore, it holds true that 

V(i,i) e r 1 x r 2 , v(t,x) e [o,r] x R fe , Vs e [t,T], Y^ mo = v^ mo {s,xt' x ). 

On the other hand, we note that Y n ' m ° is the value function of a zero-sum Dynkin game (see Appendix A4), 
i.e., it satisfies, for all s < T 

Y ir,m = ess sup CT > s ess inf T > s E[/; AT f^(r,X^, (Y^) {id)e ^x^, Z^ m °)dr+ 

max fc£(rlr ,{i;^™» -gJo-,Xt>*)}l [a<T] + {Y?>™° Awm^-^Y?™ - g^r, X^)}}l [T < a<T] 
+h^(X^ x )l [T=a=T] \T s }. 

(5.15) 

Thus, by Theorem 3.7 in Hamadene-Hassani (2005) ([15), it follows that v % ^ m ° is the unique viscosity solution 
for the following PDE with two obstacles. 

min^x) -max fce(rl) -i{u fe ^ mo (t,x) - g. k (t, x)}, max{tf(i, x) - v l ^ m ° (i, x) V mm, e(r2) - 3 (v l/ < m ° (*, x) -gji(t,x 
-df&{t,x) -b(t,x) T D x $(t,x) - iTr(aa T (t,x)DU(t,x))- 
/« (*, x, (v kl > m <> (t, x)) (M)er i xF 2 , a T (i, i).i? B tf(t, x))}] = ; 
#(T,x) = h ij (x). 

Next, let (i,x) € [0,T] x R fc and (p,q,M) E J-[v l3 ' m °}(t,x). As w ij '' m ° is a solution in a viscosity sense of the 
previous PDE with two obstacles then it holds that 

v ij ' m °(t,x)> max {v kj ' mo (t,x)-g.M,x)} (5.16) 
feg(ri)-* 

and 

max^^s) -V^ m °(t,x) Vmm le(r2r] (v a - m «(t,x) -g jt {t,x)); 

-V - b(t, x) T q - ±Tr(o-a T (t, x)M) - /«(t, x, (v M < m ° (t, x)) (M)erl xr = , a T (i, x). 9 )} > 0. ' "' " 
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But, for any constants a, b £ R we have a — (a V b) < a — b and thus a — (a V 6) > => a — fc > 0. Therefore, we 
have from (|5.17p . 

max{v^ m °(t,x) -min /e(r2rJ (^' m »(t,a;) -§ji(t,x)); 

-p - b(t, x) T q - ±Tr(aa T {t, x)M) - f^(t, x, {v kl > m ° (t, x)) (kil)er i xF 2 , a T (t, x).q)} > 0. 

Combining this inequality with f|5 . and since v^' m ° (T, x) = h l i (x) it follows that v l ^ ma is a viscosity superso- 
lution of the system 

min[z9(i,x) - max fe(rl) -i ( v k 3,™ (t, x) - g^fr a:)) ; max{i?(t,x) - min, e(r2) - 3 (v ll ' m °(t, x) -gji(t,x)); 
-d t $(t,x) -b{t,x) T D x d(t,x) - |Tr(o-<j T (t,x)D2 x ^(t,x))- 
/« (t, x, (t, x)) (fe , 0erl xr2 , a T (t, x).A^(*, x))}] = 0, 

tf(7» = /i^'(x). 

Since is arbitrary in T 1 x T 2 , the system of continuous functions (i> tJ '' m °)(t,.j)er 1 xr 2 is a supersolution of 

(O). ■ 



Consider now the set U„ lQ defined as follows. 

U = {u := (« ,3 ) (lJ ) 6r ixr= s.t. u is a subsolution of ^3TTJ> and V(i, j) e T 1 x T 2 , S*'-' < u lJ ' < #' m °}. 

Obviously, W mo is not empty since it contains (v'-'Wfjerixr 2 - Next for (i,x) e [0, T] x R fe and G T 1 x T 2 , 
let us set: 

m °v ij (t,x) =sup{u iJ (t,x), (w fci ) (/M ) e rixr2 
We now state the main result of this section. 

Theorem 5.3. The family (""f'^ijjgr'xr 2 does not depend on mo and is the unique continuous viscosity 
solution in the class Tl g of the system (2. 

Proof. We first note that for any (i, j) £ T 1 x T 2 , u* J < m °v t ^ < v % i> m ° . Since w lJ and w u,m ° are of polynomial 
growth, then ( mo v* J ')(i,j)er 1 xr 2 belongs to n g . 

The remaining of the proof is divided into two steps and mainly consists in adapting the Perron's method (see 
Crandall-Ishii-Lions, [5] Theorem 4.1, pp 23) to construct a viscosity solution to our general system of PDEs. To 
ease notation, we denote ( mo v lJ )(jj) e r 1 xr 2 by (« y )(i,j)er 1 xr 2 as no confusion is possible. 

Step 1. We first show that (w'-'Wjjer 1 xr 2 is a subsolution of (|2.1D . Indeed, it is clear that for any (t,x) £ 
[0,T] x R fc , v^(t,x) < v^(t,x) < V^' mo (t,x). This implies that < v ij <* < v^' m ° since, as pointed out 
previously, is use and v^' m ° is continuous. Therefore, for any x £ R fe , we have u y, *(T, x) = /i lJ (x), since 
v ij (T,x) = v l ^ m °(T,x) = h^(x). 

Next, fix £ T 1 x T 2 and let (« ij ')( i>5 -) e ri x r 2 be an arbitrary element of U mo . Then, for any (t, x) € [0, T) x R fe 
and any (p, q, M) £ J + v l '^*{t, x) we have 

f min[(«y> -L«((f> w '*) fc) ,))(t,x),max{({J«-* - ((v kl >*) k<l ) (t, x) , 

\ -p-&(M) T g-5Tr(aa T (t,x)Af)-fJ(^^ < 0. 

By definition we have v 1 < v kl and then v kl '* < v kl '* for any (k, I) £ V 1 x T 2 . Since the operators w = (w kl )k,i i-> 
u* 3 '* — L iJ ((w kl )k,i), w — (w kl )k,i H> u 4 -?* — U l i ((w kl )k,i) are decreasing, in view of the monotonicity property 
(H2) of the generator f*i, we have 

nin (<"•' ' - LV((v kl '*) k)l ))(t,x) ; max{(^> - {{v kl '*) k ,i))(t,x), 

-p - b{t, x) T q - ±Tr(aa T {t, x)M) - /« (t, x, [(w fe «-* (t, x)) (M)er i xr 2 , g«.* (t, x)], a T (t, x) 9 )}] < 



19 



where [(v fe ''*(t, x)) (fc .i )e ri x r 2 , v ij >*(t,x)} is obtained from (v kl '*(t, x))^;)^ 1 xr 2 ) by replacing v %3 >* (t, x) with v iJ '*(t,x). 
This means that (t,x) 6 [0,T) x M fe — > x) is a sub-solution of the following equation. 

iin[(™ - L^((t> fe; >*) M ))(t,x), max{( W - U i3 ((v kl '*) k>l ))(t,x), 

-p - b(t, x) T q - ^-Tr(aa T (f, x)M) - f 3 (t, a:, [(« fcL * (t, a:)) (k , l)e n x r 2 , w], a T (t, x)q)}} = 0. 

Next, relying on the lower semi continuity of the function 

(t,x,w,p,q,M) £ [0,T] x Rfc+i+i+fe x S fc i — ► min{(w - max fc/l (i; fc J>*(t, x) - £ ifc (*, a:))) , 
max[(w - mmi^j(v ll -*(t, x) + cjji{t, x))) , -p ~ b(t, x) T q - ^Tr(aa T (t, x)M) 
-f ij (t, x, [(v kl '* (t, x)) (fc , i)er i xr2 , w], a T {t, x)q)}} 

and using Lemma 4.2, in Crandall et al. (1992) Q5], PP-23), related to suprema of sub-solutions, combined with 
the above result, it holds that v l; > is a subsolution of the following equation: 



-p - b(t, x) T q - ±Tr(** T (t, x)M) - /« (t, x, [(v kl '* (t, x)) (M)er i xr2 , w(t, x)], a T (t, x)q)}] = 0, (5.18) 



mm[(w - LV{{v kl >*) kil )) (t, x), max{ (w - U i3 {(v kl >*) k j)) (t, x) 

-p-b{t,x) T 
v ij (T,x) = h ij (x). 

Since is arbitrary in T 1 x T 2 , (u l,? )(i,j)er 1 xr 2 is a subsolution of (12. 1[) . 

Step 2. In this step we use the so called Perron's method to show that (v iJ Wjjer 1 xr 2 is a viscosity supersolution 
of (|2~T1) . 

Indeed, we first note for any 6 T 1 x T 2 , 



since v i3 ' m ° is continuous and v 13 is ^sc. Therefore, for any x £ R fc , 

«y (T, x) — h ij (x) (5.19) 

thanks to v lJ (T,x) = h iJ (x) = v l3 ' m °(T,x). Next, assume that («*' ,7 )(i,j)er 1 xr 2 i s n °t a supersolution for (|2.1j) . 
Then, taking into account of f)5 . 19[) and Remark [1] there exists at least one pair such that v 1 ' 3 does not 

satisfy the viscosity supersolution property: this means that for some point (to, xo) G [0,T) x R fc there exists a 
triple (p,q,M) in J~ (vl r )(t ,x ) such that 

| mm[(vl 3 - L ij ((v kl )(k,i)er^xr^))(t 0l x a ) ; max{(^ J - J7*J ((i)f)(ij) er i x r 2 )) (*o, %o), , g 2Q , 

| -p-6(t ,x ) T g- 5Tr(crcr T (to,xo)M) - /' y (t , x , (wf(t , x ))( M)er i xF 2 , cr T (t , x ).<7)}] < 0. 

We now follow the same idea as in Crandall et al. (1992) ([5], pp.24). For any positive S, 7 and r, set ua j7 and 
B r as follows: 

u 5n (t,x) = v l J(t ,x ) + S+ (q,x- x ) +p(t - t ) + ^{{M ~ 2j)(x - x ), (x - x )), 

and 

B r := {(t,x) G [0,T] x R fc , s.t. |t -t \ + \x - x \ < r}. 
Using (|5.20p and continuity of all the data, choosing 5, 7 small enough we obtain 

mm[(vl 3 + 5 - L l3 {(v kl ) {ktl)er i xr 2))(t ,x ), max{(^ + 5 - U 13 ((v kl ) (k j )er i xF 2 ))(t , x ), 

-p - b(t, x) T q - iTr(aa T (to,x )(M - 2 7 )) (5.21) 
-/' y (t ,x , [(v kl (t , x ))(k,i)^(i,j) ,vl J (t , x ) + 6},a T (t ,xo)q)}} < 0. 
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Next, let us define the function T as follows. 

T(i, x) = min < us,j(t, x) — max (y* J — g. ,)(£, %), vaax{us tJ (t, x) — min (v\ +g~ji)(t, x), w(t,x)}>, 

^ k^i 1 " l^j J 

where, 

zu{t,x) = -p-b{t,x) T q - ^Tr{aa T (t,x)(M - 2 7 )) - f^{t,x, [(v^{t,x)) [kitift(itJ) ,ue^(t,x)],a T (t,x)q). 

First we note that from (|5.21l) . T(t ,xo) < 0, since ug n (t ,xo) = vl J (t ,x ) + S. On the other hand by the 
continuity of us n , Assumptions (H1)-(H2) on / lJ and finally the lower semi-continuity of u,, (fc, I) E T 1 x T 2 , we 
can easily check that the function T is use. Thus, for any e > 0, there exists 77 > such that for any (t, x) G B n 
we have T(£o,2;o) > T(t, x) — e. Since Y(io,xo) < 0, then choosing e small enough we deduce that T(t,x) < 
for any (t, x) £ -B^ with an appropriate 77. It follows that the function u$ n is a viscosity subsolution on B v of the 
following system. 

min{g(i,x) - max^, (t>^ - £.J (t, x), max{g(/j,x) - min^- + g jt )(t,x), 
-dtg(t,x) -b(t,x) J D x g{t,x) - \Tv{aa T [t,x)D 2 xx g{t,x)) 

-f ij {t, x, [{vf{t, x)) {kMh3) , x)],a T (t, x)D x g(t 7 x))}} = 0. 

Since, for any (k,l) e T 1 x T 2 , u*' < tj M '* and since satisfies the monotonicity condition (H2), U5, 7 is also a 
viscosity subsolution on B v of the system 

min{p(t,x) - max^, (v kj '* - g ik )(t,x); max{g(/j,x) - mm/^j (u a >* + ^7) (t, a;), 

-d t g{t,x) -b(t,x) T D x g{t,x) - ±Tr(aa T (<, x)D 2 xx g(t, x)) (5.22) 
-/«(*, x, [(««.* (t, x)) {k , imu) , x)], cr T (i, x)D x g(t, x))}} = 

Next, as (p,q,M) e J~{vl J {t Q ,x )) then 

v lj (t,x) > v*(t,x) > vl 3 (t Q ,x ) +p{t — t ) + (q.x - x ) 

+ i(Af(x - x ), (x - x )) + o{\t - to\) + o{\x - x | 2 ). 

2 

In view of the definition of us n and taking S = ^-7, it is easily seen that 

v %3 (t, x) > us. 7 (t, x), 

as soon as ^ < |x — xq < r and r small enough. We now take r < rj and consider the function u %3 \ 



u l3 (t,x) 



vaax.(v v (t, x), us,-y(t,x)), if (t,x)£B r , 
v 13 (t,x), otherwise. 



Then taking into account of (|5.22D and using Lemma 4.2 in Crandall et al. (1992) ([5), it follows that u 13 is 
a subsolution of (|5.18p . Next, as w' 3 > v ij and using once more the monotonicity assumption (H2) on / , 
(k,l) £ T 1 x T 2 , we get that [(v kl )(k,l)^(i,j)^' >: '] 1S a ^ so a subsolution of (|2.ip which belongs to H g . Thus, thanks 
to the Comparison Theorem 14.21 [(v fcj )(jfe i jwft 3 -),it ,J '] belongs also to U mo . Finally, in view of the definition of vl J , 
there exists a sequence (t n , x n , v rj (t n , x„))„>i that converges to {to,Xo,v % * (to,Xo)). This implies that 

lim(u y - v lJ ){t n ,x n ) = (u Sjl - vl J ){t ,x ) = vl J {t ,x ) + 5- v^(t ,x ) > 0. 

n 

This means that there are points (t n ,x n ) such that u lJ (t n , x n ) > v lJ (t n , x n ). But this contradicts the definition 
of v*i, since [(v )(fe,i)^(i,j))S y '] belongs also to U mo . Therefore, (^^(ijje^xr 2 is also a supersolution for (|2.1I) . 
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Now, in view of Corollary[2j (™°i ,! '')(jj)er 1 xr 2 is the unique continuous viscosity solution in the class U g of (|2.ip . 
Thus, using once more uniqueness, we deduce that ( mo 'y*' 7 )(i,j)er 1 xr 2 does not depend on too. M 

As above, let us denote by (v i ' :i )(i,j)er 1 xr 2 the family ( mo ^ I ' J )(i >i j)er 1 xr 2 - Here is the second main result of 
the paper. 

Theorem 5.4. For any € T 1 x V 2 , v 1 ^ = , i.e., (t>* J )(i,j)er 1 xr 2 * s continuous and is the unique viscosity 
solution in the class Tl g of the system \2.1\) . 

Proof. For any (i,j) € T 1 x T 2 and too 6 N we have 

v ij < v ij < v ij - m °. 

Taking the limit as Too —> oo we obtain w y = u y , for all G T 1 x T 2 . Finally, Theorem 15.31 completes the 

proof. ■ 

As a by-product of this result, we have the following theorem for the family (v l3 )(ij)£r 1 xr 2 - 

Theorem 5.5. The functions (V' : ')(ij)er 1 xr 2 are continuous, of polynomial growth and unique viscosity solution 
in the class H g of the following system of variational inequalities. For every 6 T 1 x T 2 , 

{max (t, x) - mm le{r2) -j(v d + g jt )(t,x), min |u« (t, x) - max te(rl) -, (v k: > - g. k )(t,x), 
-dtv? (t, x) - £«« (t, as) - /« (t, x, (« fei (t, x)) (M)er i xr2 , a T (t, x) (t, x)) } } = (5.23) 
v ij (T,x) = h ij (x). 

Proof. It is enough to consider (— u y )(i,_y)gr 1 xr 2 which, in view of Theorem 15 .41 is continuous, of polynomial 
growth and the unique viscosity solution of the following system. For all 6 T 1 x T 2 , 

!min '{t, x) — max ie ( r 2)-j {v l1 — g~j()(t,x), max (t, x) — min &e ( r i)-i (v kj +9 ik )(t,x), 
-d t vv(t,x) - Cv i *(t,x) + fi{t,x,(-v kl (t,x)) mer i x r*,-o- T (t,x)D x vV(t,x))}} =0 (5.24) 
v y '(T,x) = -h ij (x). 

Using now the result by Barles ([2], pp. 18), we obtain that (v IJ )(ij)er 1 xr 2 are continuous, of polynomial growth 
and unique viscosity solution in the class II g of system (|5.23p . which is the desired result. M 

Remark 2. We do not know whether or not we have u y = u 4 - 7 , (i,j) 6 T 1 x T 2 . 
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6 Appendix 

Al. Comparison of solutions of multi- dimensional BSDEs (|22j. Theorem 1, pp.135) 

Theorem 6.1. Let (Y,Z) (resp. (Y,Z)) be the solution of the k-dimensional BSDE associated with (/ := 
(.fi(t,u>,y,z)) i=lt k,Q (resp. (f := (fi(t,uj,y,z)) l=1 . k ,t)) where: 

(i) £ and £ are square integrable Tt -random variables ofW. k ; 

(ii) the functions f(t,co,y,z) and f(t,uj,y, z) defined on [0, T] xf2xR fc+fexd are WL k -valued, Lipschitz in (y,z) 
uniformly in (t,ui) and the process (f(t,ui,0,0))t<T (resp. (f(t,uj,0,0)) t <T) belongs to H 2 ' k ; 

(Hi) for any i — 1, . . . , k, the i — th component fi (resp. fi) of f (resp. f ) depends only on the i — th row of 
the matrix z. 

If there exists a constant C > 0, such that for any y,y € M. k , z,z€ M. kxd 

~ 4 yiifi^^'Vi +Vi>-- -,yt +Vk,z) - fi(t,u,yi, ...,y k ,z)) < 2 ^ l[ yi<0 ]\z -z\ 2 + C ivlf 

i—l,k i—l,k i—l,k 

where yf — max(yi, 0) and y^ = max(—yi, 0). Then for any i — 1, . . . , k, P-a.s., Y l <Y l . M 

A2. Systems of reflected BSDEs with one inter- connected barrier and their related systems of variational in- 
equalities (see e.g. [17] or [H]). 

Let J := {1, . . . , m} and let us consider the following functions: for i, j € J , 

fi : (t,x,y\...,y m ,z) € [0, T] x M. k+m+d f z (t, x, y\ . . . , y m , z) eM; 
g l0 : (t, x) e [0, T] x M fe h+ 9ij (t, i)eM(i/ j); 
hi-.xeR k ^ hi(x) e R. 

We now make the following assumptions. 

[Af]. For i G J , fi satisfies 
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(i) The function (t, x) H> /j(t, x, y 1 , . . . , y m , z) is continuous uniformly w.r.t. (If , z) := (j/ 1 , . . . , y m , z). 

(ii) The function fi is uniformly Lipschitz continuous with respect to (If, 2;) :— (y 1 ,..., y m , z), i.e., for some 
C* > 0, 

!/«(*,*, y\ . . . , y m , z) - fi{t,x,y\ ...,y m , z)\ < CQy 1 - y 1 ] + ■ ■ ■ + \y m - y m \ + \z - z\). 

(hi) The mapping (t,x) n> fi(t, x, 0, . . . , 0) is B([0,T] x R fe )-measurable and of polynomial growth i.e. it 
belongs to IT g ; 

(iv) Monotonicity . Mi G J, for any k G J~\ the mapping y k <E K. ^ fi(t, x,yi,..., y k -i,yk, Uk+i, ■ ■ ■ , Vm) is 
non-decreasing whenever the other components (t, x, y\, . . . , yk-i,yk+i, ■ ■ ■ , Vm) are fixed. 

[Ag]. (i) The function gij is jointly continuous in (t,x), non-negative, i.e., gij(t,x) > 0, V(i,x) G [0,T] x R fc and 
belongs to 

(ii) The no free loop property, for any (t, x) G [0, T] x R fc and for any sequence of indexes ii, . . . , ik such that 
i\ = it and card{i\, . . . , ik} = k — 1 we have 

g ili2 {t,x) + g i2i3 {t,x) + ■ ■ ■ + g lh _ llk (t,x) + g ikil (t,x) > 0, V(t,x) G [0,T] x R fe . 

As a convention we assume hereafter that gu(t,x) = for any (t,x) G [0,T] x R fe and (i,j) G T 1 x T 2 . 
[Ahl]. hi is continuous, belongs to II S and satisfies: 

Vx G R, hi(x) > max (hj(x) — gij((T,x)). 



[Ah2]. The function hi is continuous, belongs to n g and satisfies 



Vx G R, hi(x) > min (hj(x) + g i3 ((T,x)). 



Then, we have 



(6.1) 



Theorem 6.2. Assume that [Ah], [Ag] and [Ahl] are fulfilled. Then, there exist m triples of processes 

^ Y i;t,X jZ i;t,X jK i;t,x^. eJ that satis f y: V I £ J, 

Y\K l G S 2 , Z l G U 2 ' d , K l non-decreasing and K l Q = 0; 

Y* = hi(X!f) + Mr, XI'* X, vr, zi)dr + K T — K\ — jf Z*dB r , Vs<T 
Y* > m^j-^Yj - 9ij (s,X^)}, Vs < T 
{ tf(Y s * - max., {YJ gij {8,X*<*)})dKi = 0. 

Moreover there exist m deterministic functions (v l (t, x)) ie j continuous and belonging to U g such that: 

yse[t,T],Y^ x =v i (s,X t s '*). 

Finally (v l (t, x)) ie j is the unique solution, in the sub-class ofU g of continuous functions, of the following system 
of variational inequalities with inter-connected obstacles: V i G J 



min < Vi(t,x) — max (—gij(t,x) + Vj(t,x)), 

-d t Vi(t,x) - Cvfax) - fifaxyfax), . . . ,v m (t,x),a T (t,x)D x v i (t,x))} = ; (6 " 2) 
Vi(T,x) = hi(x). 
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Remark 3. In equations h6.1\) and \6.2\) . if instead we have required an upper barrier reflection, then one would 
have obtained a similar result which can be stated as follows. 



Assume that [Af], [Ag] and [Ah2] are fulfilled. Then there exist m triples of processes {(Y l ' t,x , Z r,t ' x , K r,t,x ))i^j 
that satisfy, for all i G J , 

Y\K l G S 2 , Z l G H 2 ' d , K l non- decreasing and Kg = 0; 

Y s * = hi(X*f) + jf fi(r, XP*,Y?, Y™, Z])dr - (K*. - Kj) - Z].dB r , t<s<T, 
K' ■ :"in, , ; {YJ ■ 9ij (s,X^)}, t<s<T, 
[ fi{Y* - .,• \)J + 9ij (s, X^)})dKl = 0. 

Moreover, there exist m deterministic functions {v l {t,x))i e j continuous and belong to U g such that: 

Y s i *>*=v i (s,X t s '*), s€[t,T\. 

Finally, (v l (t,x))i^j is the unique solution, in the subclass ofH g of continuous functions, of the following system 
of variational inequalities with interconnected obstacles. For all i G J 

max < Vi(t, x) — mm (gij(t, x) + Vj(t, x)), 

-dm (t, x) - Cv t (t, x) - fi (t, x, v 1 (t, x), . . . , v m (t, x), a T (t, x)D x v i (t,x))}=0: (6 ' 4) 
Vi(T, x) = hi(x). 

The proof of this result is obtained straightforward from Theorem 1 6. 6 A in considering the equations satisfied by 
((-Y\-Z\K*)) ieJ . 

A3. Linearization procedure of Lipschitz functions. Let / be a function from R 2 to R which with {x\,x 2 ) 
associates f{x\,x 2 ) which is Lipschitz in its arguments. Then, we can write 

f{xi,x 2 ) - / (1/1,3/2) = f(xi,x 2 ) - f(yi,x 2 ) + f(yi,x 2 ) - f(yi,y 2 ) 

= ^-v^ f{xi 't-l[ vl ' X2) ^ «o + **,-v W (vl, t-Z vl ' V2) ( x > - (6 - 5) 

= ai(a;i,»2,l/i)-(a;i - 2/1) + a 2 (x 2 ,yi,y 2 ).(x 2 - y 2 ) 
where, ai and a 2 are measurable functions and bounded, i.e., 

\ai(xi,x 2 ,y 1 )\ V \a 2 (x 2 ,y 1 ,y 2 )\ < C(f), i = 1,2, 

where, C{f) is the Lipschitz constant of /. Moreover, if / is non-decreasing with respect to x± (resp. x±) when 
x 2 (resp. x\) is fixed, then cij > 0, i = 1, 2. 

Linearizing / consists of writing / as 

f{xi,x 2 ) ~ / (j/i,j/a) = a 1 (xi,x 2 ,yi).(xi - y x ) + a 2 {x 2 ,yi,y 2 ).(x 2 - y 2 ). 

A4. Representation of a penalization scheme of two barriers reflected BSDE. 

For n > let (Y n ,Z n , K + > n ) be the solution of the following one barrier reflected BSDE. 

(Y n , Z n , K+' n ) G S 24 x U 2d x A+- 2 

Y? = £ + jf 9(s)ds - n jf(Y« - U s )+ds + K+' n - K t "<+ - jf Z n s dB s , t < T; (6.6) 
Y t n >L t , < t < T and J^(Y s n - L s )dK+' n = 
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where, the processes L and U belong to <S 2,1 , (g{s)) s <T G H 2 ' 1 , £ is square integrable and J^-measurable. 
Moreover, we require that L < U and Lt < £• Under these conditions, the solution (Y n , Z n , K + - n ) exists and is 
unique (see e.g. [H]). Next, for n > and t < T, set K?'~ = n f*(Y? - U s )+ds. Then, K~' n G „4+' 2 and 
Jo ( y ™ - r ™ v U s )dK~> n = 0. Therefore, the equation (Oil can be expressed as a BSDE with two reflecting 
barriers in the following manner. For all t < T, 



Thus, a result by Cvitanic and Karatzas [B] or Hamadene and Lepeltier [TB] allows to represent Y n as a value 
function of a Dynkin game, i.e., it holds true that for any t < T, 




(6.7) 



t 



ess sup CT > 4 ess inf T > t E[/ ( ffAr g(s)ds + L a l [a<r] + (F™ V f7 T )l[T<<r<T] + £,^[t=<t=t] \?t] 
ess inf r > t ess sup CT > t E[ J f °" AT g(s)ds + L a l [a<r] + (F™ V C/ t )1 [t < ct<t] + £1 [t=ct=t] |J" t ]. 



where r and cr are J-"-stopping times. 
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